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Abstract 

Radiation recoil (YORP) torques are shown to be extremely sensitive to small-scale surface topography, using nu- 
merical simulations. Starting from a set of "base objects" representative of the near-Earth object population, random 
realizations of three types of small-scale topography are added: Gaussian surface fluctuations, craters, and boulders. 
For each, the expected relative errors in the spin and obliquity components of the YORP torque caused by the observa- 
tionally unresolved small-scale topography are computed. Gaussian power, at angular scales below an observational 
limit, produces expected errors of order 100% if observations constrain the surface to a spherical harmonic order 
I < 10. For errors under 10%, the surface must be constrained to at least / = 20. A single crater with diameter roughly 
half the object's mean radius, placed at random locations, results in expected errors of several tens of percent. The 
errors scale with crater diameter D as for D > 0.3 and as for D < 0.3 mean radii. Objects that are identical 
except for the location of a single large crater can differ by factors of several in YORP torque, while being photometri- 
cally indistinguishable at the level of hundredths of a magnitude. Boulders placed randomly on identical base objects 
create torque errors roughly 3 times larger than do craters of the same diameter, with errors scaling as the square of 
the boulder diameter. A single boulder comparable to Yoshinodai on 25143 Itokawa, moved by as little as twice its 
own diameter, can alter the magnitude of the torque by factors of several, and change the sign of its spin component at 
all obliquities. Most of the total torque error produced by multiple unresolved craters is contributed by the handful of 
largest craters; but both large and small boulders contribute comparably to the total boulder-induced error. A YORP 
torque prediction derived from groundbased data can be expected to be in error by of order 100% due to unresolved 
topography. Small surface changes caused by slow spin-up or spin-down may have significant stochastic effects on the 
spin evolution of small bodies. For rotation periods between roughly 2 and 10 hours, these unpredictable changes may 
reverse the sign of the YORP torque. Objects in this spin regime may random-walk up and down in spin rate before 
the rubble-pile limit is exceeded and fissioning or loss of surface objects occurs. Similar behavior may be expected at 
rotation rates approaching the limiting values for tensile-strength dominated objects. 

Key words: ASTEROIDS, DYNAMICS, ASTEROIDS, ROTATION, ASTEROIDS, SURFACES, NEAR-EARTH 
OBJECTS 



1. Introduction 



Radiation-recoil eff^ects are now understood to be major drivers in the evolution of small bodies in the Solar 
System. The secular drift i n orbital elem ents produced by the Yarkovsky effect is a dominant mechanism in the spread 
of asteroid families ( Bo ttke et al. ^ 200 ih and in the replenishment of the near Earth object (NEO) population from the 
inner main belt (Bott ke et all 20021) . Torques produced by the YORP effect can compete with, a nd even dominate , 
collisional and tida l torques on main belt asteroids and NEOs smaller than a few km in diameter dRubincam , I2OOOI: 
Bottke et al. , 20061 and references therein). 

Recent observational efforts have achieved direct detection of the influence of the Yarkovsky and YORP mech- 
anisms. Radar observa t ions o f 6489 Golevka show non-gravitational acce lerations consisten t with Yarkovsky pre- 
dictions ( Cheslev et all 2003h. Photometric observations of 1862 Apollo ( Kaasalainen et al. . 2007 ) and 1620 Ge- 
2008ah . and photometric and radar observations of (54509) 2000 PH5 (subsequently named 



ographos dPurech et al. 
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YORP) ( ILowrv et all l2007t : iTavlor et all l2007h indicate accelerations in spin rate consistent with YORP models and 



not attributable to tidal torques. However, the YORP models contain substantial systematic uncertainties, because the 
objects' shapes and surface properties are not sufficiently constrained by observations. In the case of 1862 Apollo, 
models computed by Kaasalainen (2004) with slightly different rotation pole positions and shapes s pan a range of 
~ +25% around the measured acceleration. For 54509 YORP, the models computed bv lTavlor et alJ (120070 . usmg a 
range of shapes and surface roughnesses, span a factor of ~ 3 in predicted acceleration; and, in fact, they all systemat- 
ically overestimate the observed effect, by factors of 2 to vQ Considering the extensive observational effort devoted to 
these objects, the ambiguity in the results raises considerable concern as to whether YORP predictions derived from 
groundbased light curve and/or radar data can ever be sufficiently precise to permit quantitative tests of the theory. 

This concern is only deepened by the history of YORP calculations made for 25143 Itokawa before and after the 
arriva l of the Hayabusa spacecraft. Using pre-encounter sha p e models derived from rad ar observations ( Ostro et all 
1200 4') and photometric light curves i Kaasalainen et al.Ll2003b . IVokrouhhckv et al.l ([2004') predicted a rotational accel- 
eration of roughly half a part in 10"* per year, with an estimated uncertainty of +3 0%. Scheeres e t al. (2007) subse- 
quently recomputed the YO RP prediction, first from a s imilar pre-encounter model ( Gaskell et all I2OO6I) . obtaining a 
result consistent with that of IVokrouhhcky et al.l (120041) . and then from a sequence of shape models of progressively 
higher resolution based on in situ Hayabusa data. The first post-encounter shape model was found to imply a rota- 
tional deceleration, meaning that pre-encounter models had been unable to correctly predict even the sign of the effect. 
Subsequent models i ncorpo rating the later Hayabusa data resulted in YORP accelerations that did not converge as the 
resolution increased. ISchee res et al. (2007) conclude that, even at sub-meter resolution, the essential geometry of the 
interaction between the surface and the i ncoming and outgoing radiati on is not adequately represented. In the past year 
the situation has become even murkier. Sche eres and Gaskell (2008) show that the torque on 25143 Itokaw a is very 



sensitive to the position of the center of mass; and the predicted acceleration has still not been detected (Dure ch et al 



2008bl) . 



The foregoing stories motivate this paper, which asks whether extreme sensitivity to small-scale surface topogra- 
phy is a characteristic of the YORP effect itself, or is limited to a comparatively small number of unfortunately shaped 
objects. This question has pragmatic as well as predictive consequences. First, dynamically significant small scale 
topography may be observationally unresolvable. I will show below that, indeed, neglecting topographic features, 
such as moderate sized craters, that are indistinguishable from the ground can lead to errors in YORP predictions of 
order unity. As a result, shape models derived from groundbased data are unlikely ever to yield precise predictions 
of YORP torques. Second, dynamically significant small scale topography may, in some situations, be unpredictably 
changeable. I will demonstrate that a small displacement of a single large boulder can, in the right circumstances, 
completely reverse the sign of the YORP torque. This implies that a YORP cycle dRubincaml I2OOOI) can be inter- 
rupted in a major way by a comparatively minor event, and suggests that the spin evolution of small NEOs may be 
more stochastic than previously recognized. 

The theoretical foundations and basi c behavior of the YORP effect have been analyzed extensively by analytic, 
semi-analytic, and numerical means (e.g.,'ScheeresU2007bt| iNesvornv and Vokrouhlickvil2007ll2008clBreiter and MichalskaL 
I2OO8; Micheh and PaoHcchi, 2008; Mysen, 2008b^ These works have elucidated how asymmetries in the illumi- 
nated and radiating su rface couple to produce the now familiar behavior of the orbit-averaged torque components as a 
function of obliquity ( Vokrouhlickv and Capek , and have provided strategies for computing torques for realis- 

tic asteroid models. This paper adopts a purely numerical treatment, analyzing a large number of simulated objects to 
yield statistically meaningful results in the nonlinear regime where shadowing is important. 

The remainder of the paper is laid out as follows. In Section |2] I describe the calculation of YORP torques for 
simulated asteroids. I define a set of "base objects" on which the effect of small-scale topography will be tested, and 
describe how topography is added to the objects, in the form of Gaussian surface fluctuations, craters, or boulders. In 
Section[3] I demonstrate that small scale topography, at a level difficult or impossible to constrain from groundbased 
observations, can typically be expected to alter YORP torques by tens of percent or more. In particular, the magnitude, 
and even the sign, of the torque can hinge on the location of a single large crater or boulder. In SectionH] I discuss the 
implications of these results, both for accurate predictions of the YORP effect and for the YORP-driven evolution of 
small NEOs. 



' Scaling down the models by these factors was necessary to achieve the excellent fit in Fig. 2 of lLowrv et all j2007t) . as explained in the caption 
to that figure. 
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2. Method 



2.1. Calculation of Torques 

Radiation interacting with the surface of a body can impart a torque through the momentum it carries. It is concep- 
tually helpful to split the interaction into three parts: (1) the momentum deposited on the surface by arriving photons; 
(2) the recoil imparted to the surface by departing reflected photons; and (3) the recoil imparted by departing reradiated 
photons. In this paper I focus specifically on the third contrib ution, the reradiation torque. The fi rst contribution can 
be shown to cancel identically when averaged around an orbit ( Nesvorny and VokrouhliclcY , 20081). The second is pro 



portional to the albedo, and becomes insignificant for a sufliciently dark surface; and, as lNesvorny and Vokrouhlick^ 



(|2007) point out, if the reflected and reradiated intensities are both assumed isotropic, then the torque contributions 
are parallel. 

I compute the torques on simulated objects using the TACCH code, developed by the author and students at Ohio 
University. The object's surface is represented using a standard triangular tiling (Lagerros, 1996). For all results 
presented here, 6,394 vertices and 12,784 tiles are used, for an average angular resolution of approximately 7.5° (3 
tiles wide). For comparison, typical shape models derived from asteroid light curve observations have 1,000 to 2,000 
vertices and 2,000 to 4,000 tiles (e.g., Durech et al., 2007). Models obtained from radar observations and currently 
archived in the Planetary Data System (PDS) have typically around 4,000 vertices and 8,000 tiles. Using groundbased 
data alone, currently only 4 objects (1620 Geographos, 25143 Itokawa, 4179 Toutatis, and 52760) have been modeled 
at resolutions equal to or finer than that used in this paper. 

A horizon map (the maximum elevation (/r/, of all visible parts of the surface as a function of azimuth 77, measured 
relative to local west) is computed at the centroid of each tile when the surface is initially defined. Shadowing is 
handled by requiring that the Sun not be blocked by other parts of the surface for the tile to be illuminated. An 
ifluminated tile absorbs a flux Fabs given by 



abs 



FyUo [1 - rhem(yuo)] , (1) 



where F is the incident flux, yUo is the cosine of the angle between the incident flux and the surface norm al, and 



rhetn is the directional-hemispheric reflectance or hemispheric albedo. TACO computes rhem using eq. (42) of lHapke 



(12002b . assuming a single Henyey-G reenstein p hase function. For definiteness I adopt the mean Hapke parameters for 
S-class asteroids given bv .Helfenstein and Veverka (1989): a single-scattering albedo w - 0.23, and a phase function 
expanded to 6th order i n Legendre po lynomials with a width parameter ^ - -0.35. The correction for macroscopic 
surface roughness from iHapkd (11984^) is applied with a mean slope parameter 9 - 20°; this has no bearing on the 
absorbed flux but is relevant to the computation of lig ht curves (sect ion[3]l. These reflectance parameters are assumed 
to be constant over the surface. The opposition efl'ect ( Hapkel 2002 ) is ignored, as is the illumination of other parts of 
the surface by the reflected flux. 

I neglect the effects of thermal conductivity, so that the reradiated flux exactly balances the absorbed flux at all 
times. Each tile is assumed to radiate isotropically; if there is no blockage by the local horizon, the outgoing photons 
uniformly fill the upward-facing hemisphere, and the recoil force on the tile is given by 

= 2n f cos^ e sin Odd = , (2) 

Jo 3c 

where A is the tile area, c is the speed of light, and is the polar angle from the surface normal. In this case the 
recoil force is directed inward, normal to the surface. If, instead, part of the sky is blocked by the local horizon, some 
reradiated photons from the tile will strike other parts of the object, and the tile will intercept some photons reradiated 
by these other parts. Photons traded between different parts of the surface cannot exert a torque on the body until 
they are radiated into the clear sky. As a result, both the magnitude and direction of the recoil force will be changed. 
I compute an approximate correction for this effect, assuming that traded photons have no net effect on the energy 
absorbed by each tile. For the moment, suppose that the tile's local horizon were uniformly elevated an angle il/h 



"Thermophysical Asteroid Code, Obviously" 
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above the horizontal. Then for the same radiated flux, the intensity would be increased by a factor sec^ tj/h, and the 
magnitude of the recoil force would be given by 



/re 



In 



nc cos^ i/'/, 



■ cos sin 9d6 ■ 



2^absA 

3c 



1 + 



sin^ iff I, 
1 + sin i//h 



(3) 



The force would be increased, relative to the flat-horizon case, by a factor 1 + sin i//h/( l + sini/'/,) because the elevated 
horizon causes the outgoing radiation to be "beamed" into a narrower solid angle. In practice, of course, the elevation 
of the horizon is not uniform in azimuth. Therefore I replace (/r/, with the mean elevation, given by 



^ Ni,-l 



(4) 



where 77, is the azimuth of the /th point (out of Ni, - 50) in the horizon map. The magnitude of the recoil force is then 
taken to be amplified by a factor 1 + sin^ iphli^ + sini/'/,) over the flat-horizon case. The direction of /rec is taken to be 
antiparallel to the clear-sky normal ^sky, which is computed from a first-order Fourier expansion of the horizon map: 



A^sky = 



n.sky 
"sky' 



where 



^sky 



= A^surf - fn tan 



Ni,-\ 



TT y! ^h{r]d cos ?/,■ 
TV/, ^ 



- ew tan 



Ni,-l 



Tf y !/'/.('?/) sin 



(5) 



(6) 



Nsud is the surface normal, and and e„ are the unit vectors in the local north and west directions, respectively. 

The total instantaneous torque is then simply the cross product r x /rec, where r is the vector from the center of 
mass to the tile centroid, summed over tiles. The torque is averaged over the spin and orbital phase, assuming the 
periods are incommensurable. For computational simplicity, the orbit is taken to be circular. F or orbits with finite 
eccentricity e, the averaged torque can be simply scaled by a factor (1 - e^) (IScheeresl l2007bl) . 



2.2. Base Objects 

I define a set of base obje cts, on which to measure the ef fects of smaller scale topography, using the "Gaussian 
random sphere" formalism of iMuinonen and LagerrosI (Il998h . In this approach the distance to the surface from the 
origin is given by 



r(6, (f>) - exp 



(7) 



where (6, (p) are the usual spherical-coordinate angles and the function s(6, (p) is defined by an expansion in spherical 
harmonics: 

CO / 

5(0,0) = ^ J] s,,„ra0,^*). (8) 

/={) m=-l 

In this formulation, the mean radius of the body r = 1 . Since s{9, (p) is real, the coefficients si„, are constrained by the 
relation 

s,,^,„ = i-lT'sl. (9) 
For nonnegative m, each si,„ is a Gaussian random variable having zero expectation value and a variance given by 



Var[^R(ita)] = (l+^mo) 



In 
21+1 



Cr, 



Var[3(i,,„)] 



(10) 



In eq. ( fTOl l. S/j is the Kronecker delta symbol; C/ represents the coefficients of the Legendre expansio n of the log- 
radius covariance function, which also define the quantity /3 = (2/ Q)'^^ in eq. (|7]l- .Muinonen and LagerrosI (Il998h 
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Figure 1 : Coefficients of the Legendre expansion of the log-radius covariance function. Squares: as tabulated bv lMuinonen and Lagerroj fl99g) . 
Diamonds and dotted line: from the adopted fit, Eq. Hit . 
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Figure 2: Crater profile according to Eq. )121 . with depth to diameter ratio q = 0.2 and rim width 6 = D/20. 



estimate the C/ coefficients for / < 10 from shapes fitted to Ught curve observations of 14 objects. For this work I 
adopt an analytic fit to the results for their small-object subsample (4th column of their Table 3), given by 



Ci = 1.2 



(f + 0.26)2 



(11) 



■(^8 H- 90.0)106 ■ 

Fig.[T]shows the coefficients from Muinonen and Lagerrosl ( |l998|) and the fit. This fit will be used below to extrapolate 
the covariance function out as far as / = 40. For Z » 1 , Eqs. (fTOl i and (fTTT i imply C/ ~ l^'^ '^^, so that the RMS amplitudes 
Sim decline as T^'^"*. 

Each base object is a Gaussian random sphere with the expansion in Eq. ((S]) truncated at a maximum order /base- 
I adopt values of Zbase = 4,5,7,9, 12, 15,20, and, for each value, compute 6 different base objects with independent 
random realizations of the 5/„, coefficients. 

2.3. Addition of Small- Scale Topography 

The base objects are intended to represent asteroids whose shapes are constrained by observations down to some 
minimum angular scale. Different types of small-scale topographic structure can be added to the base objects. In 
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Figure 3: Typical variation in YORP torques caused by Gaussian topograpliic power added at smaller scales to the base objects, (a) Base object 
with /(jase = 9, topography added to /,^iix = 18. Top- and bottom-most curves coiTespond to the objects shown in Fig. |4] (b) Base object with 
'base = 20, topography added to /^ax = 40. Top panels show the spin component, and bottom panels the obHquity component, of the torque, 
normalized by the moment of inertia about the short (spin) axis, plotted against obliquity. Individual curves correspond to ten different random 
realizations of the added topography. 



Sectional measure the separate effects on YORP of each of the following: 

1 . Gaussian power at smaller scales. Here the expansion in Eq. ([8]l is simply extended from /base to some higher 
maximum order /max- The si,„ coefficients are identical to those for the base object for I < /ba,se, and chosen 
randomly according to Eq. ( fTOl) for /base < I ^ ^max- 

Craters. One or more circular craters are placed randomly on the surface. The crater depth z as a function of 
radius r from its center is given by 



2. 



z(r) = qD 



l,if r <D/2; 



exp 



ifD/2 <r<D/2 + 3d. 



(12) 



In Eq. (fT2]) . D is the crater diameter, q is the depth-to-diameter ratio, and 6 is the width of the raised rim. I 
adopt q = 0.2 and 6 = D/20, based on NEAR- Shoemak er observations of 253 Mathilde (Tho mas et aL , 
and 433 Eros ( Wverka et all 2001 : Thomas et al. . 2002) and laboratory impact experiments (iNakamura . 



1999) 



2002 



Housen and Holsapplei 2003). The resulting profile is shown in Fig. |2] The vertical displacement z is in the 



direction of the mean surface normal of the area the crater will occupy. Prexisting topography is first flattened to 
the level of the mean surface and then overwritten by the crater Multiple craters are added sequentially so that 
overlapping features are reasonably realistic. 
3. Boulders. One or more boulders can also be placed on the surface. As a crude representation of a boulder, the tile 
vertices falling within a given radius of the boulder center are raised, creating an approximately round, mesa-like 
feature with jagged edges. As with craters, the existing topography is first flattened so that the boulder's upper 
surface is featureless. 

For each object, the center of mass and inertia tensor are calculated after the small-scale topography is added, assuming 
a homogeneous density p. The coordinate system is shifted and rotated to align the origin with the center of mass 
and the coordinate axes with the principal axes of the body. All objects are assumed to be rotating about the shortest 
principal axis. 



3. Results 

3.1. Gaussian Power at Smaller Scales 

Any set of observations will be able to constrain the topographic power down to some minimum length scale or 
maximum harmonic order. As a straightforward indicator of the effect of neglected higher orders, I add power to each 
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Figure 4: Identical views of the two Zbasc = 9 objects from Fig. [3^ with the (left) smallest and (right) largest torques. Line of sight is in the plane 
containing the short and middle axes of the body, 45° from each, and the illumination is from the right at 45° phase angle. 



base object to a maximum order Zmax = 2Zbase, and measure the error in torque incurred by ignoring topography down 
to half the observationally constrained scale. 

Fig.[3]shows typical results for two of the 42 base objects, one with Zbase = 9 and one with /base = 20, corresponding 
to minimum angular scales of 40° and 18°, respectively. The curves show the averaged torque components affecting 
the spin rate (r ,) and the obliquity (t^), normalized by the principal moment of inertia about the short axis plotted 
against obliquityH Different curves indicate the 10 random realizations of the added topography. Clearly, knowing 
the topography only to order /base - 9 (Fig. [3^) is not sufficient to determine the torque to even order-unity accuracy. 
The range in r, spans nearly a fact or of 3 at obliquity e = 0°, and more than a factor of 5 at e = 90°. The location 
of the "Slivan states" dSlivanl 2002 ). which require = 0, varies over 15°; and the obliquity torque in the vicinity 
of the Slivan states varies over a factor of 5. In contrast, the /base - 20 objects (Fig.IlJ)), to which power has been 
added to /^ax = 40, have largely concordant torques. Evidently, at sufficiently small angular scales, the neglect of 
Gaussian power at still smaller scales is not too serious, as long as it follows a simple extrapolation of the covariance 
coefficients. 

As a measure of how wrong a YORP prediction is likely to be due to neglect of the smaller scale topography, I 
compute, at each value of obliquity, the standard deviation in each torque component over the 10 trials, then average 
over obliquity; this quantity is normalized to the mean absolute value of that component, also averaged over obliquity. 
The result is an obliquity-independent measure of the expected relative error For the /base - 9 case in Fig. [3^, the 
expected relative errors are 38% in Tj and 32% in t^. Keep in mind that these are Icr errors, implying 95% confidence 
bands (for a Gaussian distribution) with full widths of 152% and 128%, consistent with the order-unity variation seen 
in the figure. For the /base = 20 case, the relative Icr errors are 7% and 3% in r, and respectively. I will use this 
same approach to quantify the effects of craters and boulders below. 

The most remarkable thing about the results in Fig. [3^ is that the individual objects, despite differing by factors 
of a few in their YORP torques, appear nearly identical. Fig. |4] shows a representative view of the two objects in 
Fig. [3^ having the most discordant torques. To assess whether these objects would be distinguishable photometrically, 
I compute optical light curves over a grid of 40 rotation pole directions and 10 phase angles from 0° to 90°. I find that 
the RMS photometric difference between these two objects is < 0.05 magnitudes at all phase angles. The maximum 
difference is < 0.1 mag at phase angles < 50°, and is > 0.2 mag only at phase angles of 75° or higher, making it 
unlikely that these objects would be distinguished by groundbased observations. 

In Fig. |5] I show the relative 1-sigma errors for all of the base objects. The scatter at each /base reflects the base 
objects' different YORP susceptibilities. One can easily see that order-unity errors in the torque should be expected for 
'base < 10. To be reasonably confident of errors under 10%, the surface must be known to at least /base = 20. Owing to 



'Numerical results for normalized torques are given in units of 1.52 X 10 '^(p/2gcm '(r/1 km) ^(a/1 AU) ^(1 - e^) ''^ s ^, where p and 
r are the constant density and mean radius of the body, and a and e are the orbital semimajor axis and eccentricity. 
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Figure 5: Expected relative Icr error in the YORP torque caused by neglecting Gaussian topographic power at /base < / < 2/basc- Top: eiTor in the 
spin component; bottom: error in the obliquity component. The relative error is the standard deviation in each component over 10 random trials, 
averaged over obliquity, and normalized to the mean absolute value of that component. Multiple symbols correspond to 6 different base objects at 
each /base- 



the steepness of the extrapolated covariance function, most of the change in the torque is actually produced by orders 
close to /base- Fig.|6]shows the relative errors for the Zba,ie - 12 objects with power added to /max = 13, 14, 16, 18,24, 
and 36. There is no discernible increase in the error for /max > 18, and the single order 1-13 seems to be responsible 
for roughly half of the asymptotic value. 



3.2. Craters 

Actual small-scale topography is likely to be non-Gaussian. Impact craters are one example of Poisson-distributed 
features not well described by a truncated spherical-ha rmonic expans ion. The c rater size dis tribution has b een esti- 
mated for 951 Gaspra (Chapman et al., 1996b; Green berg et all 1 1994, 243 Ida dBelton et al.L 1994; Chapm an et al 



1996ah . 253 Mathilde (Chapm an et all 1 19991) . and 433 Eros dChapman et all l2002h . To lowest order, the cumula 



tive distributio n of craters w i th dia meters larger than D, N{D), is reasonably well approximated by a power law, 



N(D) ~ D-^ dO'Brien et all l2006h . Though there are significant deviations in the logarithmic slope for different 



objects, the known distributions all describe a steeply diminishin g number of large c raters. The largest craters have 



diameters ranging from 0.35 to 1.2 times the bodies' mean radii ( Thomas et allll999h . 



I return to the original base objects, removing the small-scale Gaussian power added in Section 13.11 in order 
to measure separately the effect of craters on the YORP torque components. I start by randomly placing a single, 
modest-sized crater of diameter D - 0.3 (in units of the body mean radius) on the base objects. Fig. |7] shows the 
resulting relative Icr errors, calculated as in Sec. 13.11 with ten random crater placements for each object. A single 
crater can clearly alter the spin and obliquity torques by tens of percent, even though the crater itself, including the 
raised rim, occupies only ~ 0.5% of the total surface. The effect increases with /base at values /base ^ 9. A similar, 
but less pronounced, increase is also seen in the Gaussian case (Fig.|5]l. The single crater actually produces less error 
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Figure 6: Expected relative Icr error, as in Fig. [5] for the /(,asc = 12 base objects with Gaussian power added to /,nax = 13, 14, 16, 18, 24, 36. There 
is no strong dependence on /^ax for 'max > 18, indicating that orders just above the cutoff /base most important. 



for /base ^ 9 than does small-scale Gaussian power added to the same base object. But the effect of the crater for 
'base > 10 is larger than the Gaussian case and, to within the statistical errors, independent of the base order This may 
be due either to the declining ampUtudes C; at higher I or to the fact that I - 9 corresponds roughly to the scale of the 
crater itself. Evidently the crater couples to other mesoscale surface features of roughly its own size more effectively 
than to the large-scale structure at smaller /. 

Fig. shows the normalized torque as a function of obliquity for ten random placements of single craters on 
the same /base = 20 base object used in Fig. |3]3. The effect of the crater is about twice that of the added Gaussian 
topography: the relative errors in and are 14% and 6%, respectively, which also represent median results over all 
base objects. Fig.jS]? shows the effect of doubling the crater size to D = 0.6 on the same base object. The relative errors 
increase to 56% and 14%. Comparing the extreme examples (the highest and lowest curves), the spin components 
differ by a factor of 4.4 at e = 90° and actually have opposite signs at e = 0°. Fig.|9]shows these two most discordant 
objects, which turn out to be cratered at almost diametrically opposite positions. 

Cratering a base object moves the center of mass and slightly reorients the principal axes. To determine whether 
these adjustments contribute to the change in torque, I recompute the torques for 60 of the cratered and bouldered 
(see Section [33T l objects, without shifting to the new center of mass or reorienting to the new principal axes. For 
either craters or boulders, the results differ, on average, by only a few percent, indicating that it is really the surface 
topography that is responsible for the large expected torque errors. 

The objects, even those with single large (D - 0.6) craters, are so photometrically similar that it is unlikely that 
they would be distinguishable in current light curve programs. To assess this possibility, I take all 60 of the /base = 20 
cratered objects and compute light curves over the same grid of 40 orientations and 10 phase angles used in section 
13.11 I then calculate the probability that two objects, differing only in crater location, would be observed at a given 
moment to have a magnitude difference > A. The result, as a function of phase angle, is shown in Fig. \TQ\ Note 
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Figure 7: Expected relative Icr standard deviation in the YORP torque, as in Fig.|5] caused by a single randomly placed crater with diameter 
D = 0.3 added to different base objects. 



that, unless substantial data are obtained at phase angles > 70°, one expects that no more than 1% of the total data 
will differ by more than 0.05 mag, and no more than 0.01% by more than 0.1 mag. In practice, observations at phase 
angles above 45° are rare in light curve programs. 

Fig.[TT]shows the scaling of the relative torque errors with crater size. The diamonds and plus signs show results 
for single-cratered objects with /base =12 and 20, respectively. One might naively expect the torque error to scale with 
the crater area (~ D^)\ and Fig.[TT]shows that this is not a bad description for D > 0.3. However, at smaller diameters 
the scaling is much steeper, with the relative error roughly oc D^. This nonlinear response may be attributable to 
shadowing, although the same steepening does not occur for boulders (Section FS.SI below). The squares and crosses 
in Fig.[TT](slightly offset for clarity) indicate the results for objects with four D = 0.3 craters each. Since the changes 
in torque produced by each added crater are random and uncorrelated, one should expect the total change to scale as 
the square root of the number of craters. This expectation is confirmed by the factor 2 increase in the relative errors 
over the single-crater case. The implications of these scaling laws for the expected accuracy of YORP predictions on 
realsitically cratered objects are discussed in section|4]below. 



3.3. Boulders 

Boulders are a second example of Poisson-distributed topography. The surface of 25143 Itokawa has been found 
to be dominated by boulders, with a cumulative size distribution well represented by N( D) oc ov er a factor 10 in 
diam eter; the largest boulde r, Yoshinodai, is fully one-tenth the size of its parent body (iSaito et al lllOOej). For 433 
Eros, IChapman et alj (120021) obtain a result consistent with N(D) oc D These size distributions are substantially 
steeper than those for craters; thus one can expect the cumulative effect of moderate sized boulders to be more 
relatively important than that of moderate sized craters. 

I proceed with the same strategy to assess the influence of boulders on YORP torques. I add a single mesa-shaped 
boulder, with diameter D = 0.3 (in units of the body mean radius) and height h = D/2, at random locations to the base 
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Figure 8: Typical variation in YORP torques caused by single craters of diameter (a) D = 0.3, (b) D = 0.6, randomly placed on the same base 
object used in Fig.[3jj. Top- and bottom-most curves con'espond to the objects shown in Fig.|9] 




objects. The base objects have mean longest diameters of approximately 3.1, so the size of these boulders relative 
to their parent bodies approximately match es that of Yoshinoda i. The height-to-diameter ratio corresponds to typical 



values measured for boulders on 433 Eros (iThomas et alil2002h . 



Fig. [12] shows the expected relative error in the torque components caused by the single large boulder The 
effect is quite dramatic: typical Icr errors are 10% to 100% in the spin component, and 5% to 50% in the obliquity 
component. This is roughly a factor 3 larger than the effect of craters of the same diameter, though the difference may 
be attributable to the simple fact that boulders can be proportionally taller than craters are deep (here, by a factor of 
2.5). As in the case of craters, there is a weak increase in the mean error for values of Zbase < 10, but little difference 
is seen at higher order The scaling of the errors with boulder size is shown in Fig. [13] The boulders are all similarly 
shaped, with h - D/2. The relative error scales quite cleanly as D^, i.e., proportional to the boulder surface area. 

Examples of the torque components as a function of obliquity are shown in Fig. [T4h . for 10 random D = 0.3 
boulder placements on the same Zbase = 20 object used in Figs. [3J) and[8] This case again represents a typical amount 
of topographically induced variation. Boulders have a tendency to shift the curves for Tj / /, up and down while altering 
the overall amplitude of those for t^/Is- Note that the sign of the spin torque can be flipped by the location of a single 
boulder, even for this unexceptional object. A more extreme example is seen in Fig. [T4b . which shows the same 
curves for one of the other five /base - 20 objects. Here the sign of the spin torque at all obliquities can be reversed by 
the location of a single Yoshinodai-sized boulder The two most discordant objects are shown in Fig. [15] The boulder 
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Figure 10: Probability contours for observing a photometric difference A or larger, as a function of phase, between pairs of objects that dilfer only 
in the location of a single crater of diameter D = 0.6. Results shown are averaged over the 6 base objects with /base = 20. 



location differs by barely twice its own diameter; yet the object on the left will inexorably spin up while the object on 
the right will spin down. 



4. Discussion 

The above results indicate that sensitivity to small-scale topography is an inherent characteristic of the YORP 
effect. Since YORP operates by virtue of the non-cancellation of many independent torque components that vary 
across the surface, over the rotation of the body, and around the orbit, this conclusion is perhaps not surprising; 
however, it does impose strong limits on the predictability of YORP torques, and has important consequences for the 
spin evolution of small bodies. 

If the surface is describable by a spherical harmonic expansion whose coefficients are Gaussian-distributed, the 
topography must be known accurately at least to harmonic order I - 10 to allow a YORP prediction of even the correct 
sign. To achieve 10% accuracy requires the surface to be kn own at least to / - 20, if t he topography diminishes with 
/ in a manner consistent with a simple extrapolation of the iMuinonen and LagerrosI (fl998) covariance coefficients 
(Eq. ifTTI ). Few objects can be modeled reliably at this resolution. Furthermore, actual small scale topography is not 
Gaussian and not small in amplitude. The above results show that Poisson-distributed features such as craters and 
boulders, which produce substantial shadowing and which are not well revealed by groundbased observations, can 
have significant effects. A single large feature can alter YORP torques by tens of percent or more. 

1 found above that the relative torque error scales as cr ~ D^, where [3^2 for boulders and 2 < /3 <3 for craters. 
Using these numerical results and the observed size distributions, one can estimate the relative influence of large and 
small craters and boulders. The observed cumulative counts are reasonably well represented by N(D) ~ D ", with 
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Figure 1 1 : Expected relative en'or in torque, as in Fig.0 caused by single craters randomly placed on base objects with /base = 12 {plus signs) and 
Wse = 20 (diamonds), as a function of crater diameter D. Crosses and squares (offset horizontally) show results for 4 randomly placed craters each 
of diameter D = 0.3, on objects with /base = 12, 20 respectively. 



3 < a < 4 for boulders and I < a < 2 for craters0 Assuming the effects of individual boulders or craters to be 
uncorrected, one expects in either case that the variance in torque produced by features with diameters between Dmin 
and Dmax will scale according to 



r 



dN_ 
dD' 



D^'^dD oc D 



ip-a 
max 



D 



{a + ip). 



(13) 



For cratered objects, Eq. ( fT3l l implies that at least 75% of the total variance (87% of the expected cr) is produced 
by the 4 biggest craters. Fig.[TT]shows that the expected error from the random location of the single largest crater is 
already approaching 100% at D^ax = 0.6, close to the observed maxi mum diameter on Ero s but rather modest com- 
pared to the values of 0.89 and 1.2 for Ida and Mathilde, respectively ( Thomas et al. . 19991) . These very large craters 
could possibly give rise to noticeable features in high phase-angle light curves; however, l ight curve inversion per- 
forms poorly in recovering surf ace concavities because the inversion problem is non-unique ( Kaasalainen et all 1992 ; 
Kaasalainen and Torppal l200lh . The convex, pre-encounter model of the Rosetta target 2867 Steins dLamv et al. , 
2008 ) resembles the in situ images jf| but by construction misses several large craters, including one with an apparent 
diameter D «; 0.9. A non-convex inversion probably would not have revealed this crater, as the light curve data are 
limited to phase angles < 42°. Convex models, in general, may be useful for determining overall surface reflectance 
and thermal properties, and may give accurate predictions for Yarkovsky accelerations. But a YORP prediction based 
on a convex model is likely to be wrong by of order 100%. 



*As before, diameters are implicitly in units of the body mean radius, 
^http: //www. esa.int/SPECIALS/Rosetta/SEMN MY04KKF_0.html 
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Figure 12: Expected relative Icr standard deviation in the YORP torque, as in Fig.|5] caused by a single randomly placed boulder with diameter 
D = 0.3 added to different base objects. 



Boulders, because of their steeper size distribution, behave differently: at most half of the total variance is pro- 
duced by boulders larger than half the maximum size. In fact, if the distribution is as steep as on Eros (a - 4), then 
the variance is formally divergent as D — » 0. Of course, distributions dN/dD ~ D ' or steeper are unphysical in 
this limit; but the scaling indicates that the YORP torque can be affected as much by the surface distribution of small 
boulders as by that of large ones. The results of section[3]also show that single boulders may have significantly greater 
effect on the torque than craters of the same size; a single Yoshinodai-sized object can easily alter the spin torque by 
tens of percent, and in some cases change its sign. The clear implication is that the surface boulder distribution can be 
critical in determining the actual consequences of YORP. 

Thus, one should expect uncertainties of order unity in predictions of the YORP effect made from models derived 
from remote observations, just from the unconstrained nature of the small-scale surface topography. The true situation 
is actually even worse, since 1 have neglected other issues, such as variations in a l bedo, surface roughness, and internal 
density, which probably have comparable effects. Indeed. IScheeres and Gaskelll (l2008r) find that shifting the center of 
mass of Itokawa by only 15 m can reverse the sign of the spin torque. At present, there have been 4 attempts to predict 
the YORP acceleration from groundbased data and compare with direct measurements. In one case (1682 Apollo), 
the correct magnitude and sign were obtained, with an uncertainty of ±25%; in a second (54509 YORP), models 
gave the correct sign but the wrong magnitude by a factor of 2 to 7. In the third (25143 Itokawa), pre-Hayabusa 
models evidently obtained the wrong sign (compared to later models), but the predicted acceleration has still not been 
detected. And in the fourth (1620 Geographos), torques computed for different convex shape models bracketed the 
observed value, but varied over a factor of 5. This situation is about what we should expect. Generally speaking, a 
YORP prediction computed from a groundbased model will be accurate only to within factors of order unity, and will 
have no better than an 80% chance of having the correct sign. 

The sensitivity of reradiation torques to boulders also implies that there may be substantial stochasticity in the spin 
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Figure 13: Expected relative en'or in torque, as in Fig. 1 121 caused by single boulders, as a function of boulder diameter D. Symbols have the same 
meaning as in Fig. II II 



evolution of small bodies under the influence of YORP. In the standard picture of the YORP cycle d Rubincam , 2000l) . 
an object can gradually spin up, simultaneously evolving in obliquity, as determined by the now familiar "YORP 
curves" (e.g., Figs. [3] E] and[T4]i. In most cases, spin-up continues only until reaching an orientation regime where 
Ts < 0. The object then spins down, asymptotically approaching the stable fixed point in obliquity. Objects with 
substantial power in the higher YORP orders (Nesvorny and Vokrouhlicky, 2007, 2008) can sometimes have stable 
fixed points at obliquities where Tj > 0; examples of this behavior are seen in Figs.[8]3 andfT4b. along with cases where 
Tj > at all obliquities. These objects would presumably spin up indefinitely. But a new possibility raised here is 
that a seemingly minor topographic change, such as the movement of a boulder, can stochastically shift the evolution 
onto a new set of YORP curves. One would expect such events to occur well before catastrophic fissioning or even 
unbinding of surface material, as the loca l effective gravity vectors change with the accelerating spin. Calculations 
for contact-binary configurations (iScheere s. 2007a) and numerical experiments on rubble piles ( Walsh et al, 200i 
Harris et al. . l2009h confirm that surface changes occur at rotation rates slower, by factors of several, than the nominal 



centrifugal limit for self-gravitating spheres. The latter corresponds to the observed minimum rotation periods of large 
asteroids, at approximately 2 hours. Thus a broad range of spin periods, from 2 to perhaps 10 hours, may correspond 
to a "scattering zone", in which small topographical changes stochastically reverse the YORP torques. An object 
being spun up by YORP may scatter multiple times through motion or loss of surface boulders, random-walking up 
and down in spin rate, before finally crossing the actual spin limit and either fissioning or shedding all loose objects. 
It is extremely in teresting that the observed distribution of NEOs in the diameter-rotation period plane (e.g., Fig. 1 of 



Holsapplei l2007l Statler et al. 2009, in preparation) shows the sharp cutoff at 2 hours becoming notably indistinct at 



diameters < 1 km, where YORP is expected to b ecome dominant. Similarly, small, fast-rotating objects approaching 
the higher spin limit imposed by tensile strength dHolsapplelEoOTl) may wander up and down in spin rate, driven not 
by catastrophic fragmentation, but merely by a subtle rearrangement of the surface that reverses the YORP eff'ect. 
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Figure 14: Typical variation in YORP torques caused by single boulders of diameter D = 0.3, randomly placed on {a) the same base object used 
in Figs.[3j) and[8] and (b) a second object with /b^se = 20, showing more extreme variation. Top- and bottom-most curves correspond to the objects 
shown in Fig. [15] Note in (b) the possibility of reversing the sign of the spin torque at all obliquities simply by repositioning one boulder. 




Figure 15: Identical views of the two objects from Fig. 114b with the {left) most positive and (right) most negative spin torques. Line of sight is in 
the plane containing the short and middle axes of the body, 45° from each, and the illumination is from the right at 30° phase angle. Moving the 
boulder by only twice its own diameter reverses the sign of the spin torque. 
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